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ABSTRACT 

O . 

O ■ We present an analytical and phenomenological model for metal enrichment 

in halos based on hierarchical structure formation. This model assumes that 



astration of normal stellar populations along with Type II supernovae (SNe II) 
already occurs at very high redshift. It focuses on the regime of [Fe/H] < —1 
where SNe la are not major contributors. The regime of [Fe/H] > — 1 where 
CN \ SNe la are major contributors is also discussed in general. For halos that are 

! not disrupted by SN II explosions, the chemical evolution of the gas and stars is 

explicitly determined by the rate of gas infall as compared with the astration rate 
\ and the corresponding rate of metal production by SNe II per H atom in the gas. 

^ | This model provides a good description of the data on [Fe/H] for damped Ly« 

systems over a wide range of redshift 0.5 < z < 5. For all halos not disrupted by 
SN II explosions, if there is a cessation of gas infall, the metallicities of stars follow 
a bimodal distribution. This distribution is characterized by a sharp peak at the 
value of [Fe/H] corresponding to the time of infall cessation and by a broad peak 
at a higher value of [Fe/H] corresponding to the subsequent period of astration 
during which the bulk of the remaining gas forms stars. Such a distribution may 
be compared to that observed for the Galactic halo stars. If the gas in a halo 
. ^ | is rapidly lost upon cessation of infall, then an assemblage of stars with a very 

sharply-defined [Fe/H] value will be left behind. This assemblage of stars may be 
accreted by a larger system and become a globular cluster of the larger system. 
We also discuss the masses and metallicities of the globular clusters in this model 
considering the criterion for the onset of astration and possible disruption of low- 
mass halos by SN II explosions. The implications of possible globular clusters 
with very low metallicities of [Fe/H] < —3 are discussed. 

Subject headings: galaxies: abundances — galaxies: evolution — intergalactic 
medium 
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1. 



Introduction 



We present an analytical and phenomenological model for chemical evolution of galaxies 
based on hierarchical structure formation. This model provides a possible explanation for the 
metal abundances observed in damped Lja (DLA) systems and the metallicity distribution 
of stars in galaxies. It also suggests a possible mechanism for the formation of globular 
clusters. Measurement of "metallicities" of DLA systems by Lu et al. (1996) and further 
studies by Prochaska & Wolfe (2000, 2002) showed that these systems have a large dispersion 
in [Fe/H] = log (Fe/H) — log (Fe/H) Q for a given redshift z and that there appears to be a 
baseline enrichment at the level of [Fe/H] ~ —3 over a wide range of redshift 1.5 < z < 4.5. 
The work by Prochaska et al. (2003) provided a much more extensive data base over a 
wider range of redshift 0.5 < z < 5. All these data exhibited the same behavior. We earlier 
found that the data on DLA systems could be explained by a simple closed-system model 
(Wasserburg & Qian 2000b and Fig. 1 here), where the number abundance (Fe/H) of Fe 
relative to H atoms in a DLA system is determined by 



In the above equation, (Fe/H)p ~ lO _3 (Fe/H) ([Fe/H]p ~ —3) is the initial Fe abundance 
corresponding to a postulated prompt inventory of metals in the intergalactic medium (IGM), 
Ap e is the rate of Fe production by Type II supernovae (SNe II) per H atom in the gas of 
the DLA system and is taken to be ~ (30 Gyr) -1 based on estimates for the Galaxy (see 
§2), t(z) is the age of the universe at the redshift z of the system, and t* is the "turn-on" 
time corresponding to the onset of astration and Fe production in the system. It was found 
that the available data on [Fe/H] for essentially all the DLA systems lie between the prompt 
inventory value and the upper bound obtained for t* = 0. The large dispersion in [Fe/H] at 
a given z could be explained by different t* for different DLA systems. It was found that for 
the above closed-system model, most systems at a given z have t* in a range close to t(z) 
(i.e., they "turned on" long after the big bang so that their baryonic matter had to be stored 
in the IGM for this long period). 

The astrophysical implications of the late turn-on of most DLA systems in the above 
closed-system model were explored by us in a recent study based on the cold dark matter 
(CDM) model of hierarchical structure formation (Qian & Wasserburg 2003b). We there 
associated the evolution of a DLA system with the growth of a halo through infall of CDM 
and gas. The late turn-on was reinterpreted in terms of two factors. First, a halo must 
achieve some threshold mass before star formation could occur. As a halo associated with a 
low-cr density fluctuation reaches a given mass at a later time, Fe production in such a halo 
would be delayed simply due to the later onset of star formation. Second, subsequent to 
the onset of star formation, the growth of (Fe/H) in the gas of a halo is determined by the 



(Fe/H) = (Fe/H) p + A Fe [t(^) 




(1) 
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competition between production by SNe II and dilution by gas infall. The uniform growth 
with time does not set in until infall ceases. Thus, the late "turn-on" does not require the 
introduction of an arbitrary delay time t* but is a direct result of the infall rates in the CDM 
model of hierarchical structure formation. 

We here extend our model for chemical evolution of DLA systems to a more general 
discussion of galactic chemical evolution. In §2 we review the consequences of astration 
concurrent with infall for the chemical evolution of gas in a halo that is not disrupted by SNe 
II. We discuss the effects of infall cessation, particularly the resultant metallicity distribution 
of stars in §3. The connection with the CDM model of hierarchical structure formation is 
described in §4. A comparison of our chemical evolution model with the data on DLA 
systems is presented in §5. Possible scenarios of globular cluster formation in connection 
with this model are discussed and conclusions given in §6. 

2. Chemical Evolution of Gas in an Accreting Halo without Disruption 

We first consider the chemical evolution of gas in a halo that grows according to hi- 
erarchical structure formation. We assume that CDM and gas are fed into the halo at a 
fixed mass ratio through infall. For simplicity, the infalling gas from the IGM is assumed to 
be pure big bang debris (i.e., metal-free). Observations show that the IGM is significantly 
enriched with "metallicities" corresponding to [Fe/H] ~ —3 over a wide range of redshift 
1.5 < z < 5.5 (e.g., Songaila 2001; Pettini et al. 2003; see also Simcoe, Sargent, & Rauch 
2002, 2004; Schaye et al. 2003; Aguirre et al. 2004). Thus, the assumption of metal-free in- 
falling gas is not necessarily valid. However, possible metallicities at the level of [Fe/H] ~ —3 
for the infalling gas are soon overwhelmed by enrichments associated with astration in the 
halo. Insofar as the evolution of the halo substantially enhances the metal content of its gas 
to well above [Fe/H] ~ —3, neglecting the metallicity of the infalling gas is justified. 

The history of star formation in the halo is crucial to the chemical evolution of the 
gas. The issue of what the first generation of stars were like has long been a matter of 
intense interest. Numerical simulations of collapse and fragmentation of metal-free gas clouds 
indicate that the first stars (or baryonic aggregates) were very massive. Formation of normal 
stars may not have occurred until a "metallicity" corresponding to [Fe/H] ~ —4 was reached 
in the IGM (see Bromm & Larson 2004, and references therein). We do not consider the 
effect of this threshold metallicity that may possibly be required for normal astration. We 
assume that normal astration proceeds at a constant rate per H atom in the gas once the 
halo reaches some critical mass (see §4). From the discovery of a low-mass Galactic halo 
star with [Fe/H] = —5.3 (Christlieb et al. 2002), it is evident that some low-mass stars were 
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formed when the level of heavy metals in the IGM was extremely low, far below the threshold 
metallicity for normal astration suggested by theoretical models (Bromm et al. 2001). 

It has long been recognized that stars in the Galactic halo form a distinct population. 
Extensive observational work has established that these stars have distinctively low [Fe/H] 
ranging from typically ~ —2 (e.g, Carney et al. 1996) down to —5.3 (see Fig. 5 in Christlieb 
2003). There is a clear break at [Fe/H] ~ —1 observed in the Galactic evolution of [Ba/Fe] 
and [O/Fe] with [Fe/H]. For [Ba/Fe], this break is due to the lack of major s-process contri- 
butions from asymptotic giant branch stars to the Ba in the interstellar medium (ISM) for 
[Fe/H] < —1 (cf., Busso, Gallino, & Wasserburg 1999), while for [O/Fe], it is due to the lack 
of major Fe contributions from SNe la for [Fe/H] < —1 (Tinsley 1979; see also Kobayashi et 
al. 1998). In addition, the abundance patterns of the elements from the Fe group down to 
O observed in Galactic halo stars are characteristic of the nucleosynthetic yields of SNe II. 
Thus, the overall abundances of "metals" in the early Galaxy corresponding to [Fe/H] < — 1 
in the ISM were clearly dominated by contributions from SNe II with massive (~ 10-100 M ) 
progenitors. We note that the marked change at [Fe/H] ~ —3 observed for the abundances 
of heavy r-process nuclei in Galactic halo stars (e.g., McWilliam et al. 1995; Sneden et 
al. 1996; Westin et al. 2000; Burris et al. 2000) appears to require a source other than 
SNe II for these nuclei (Qian & Wasserburg 2003a). However, this r-process source cannot 
produce any significant amount of the elements from the Fe group down to O (e.g., Qian & 
Wasserburg 2002), and therefore, does not affect our discussion of the chemical evolution of 
these elements here. 

Based on the above discussion, we only consider metal production by SNe II subsequent 
to the onset of astration during the early evolution of a halo. The general effects of SNe 
la that contribute mainly the Fe group elements at [Fe/H] > —1 will be discussed in §3. 
For convenience, we use Fe as a representative metal element. At [Fe/H] < —1, the other 
elements (from the Fe group down to O) essentially evolve in proportion to Fe as specified 
by the effective yield pattern of SNe II (see Table 1 in Qian & Wasserburg 2002). If gas 
in the halo is homogeneous and is not expelled by explosions of SNe II, the basic equation 
governing its chemical evolution (Qian & Wasserburg 2003b) may be written as: 

<i(Fe/H) P Fc 1 d(R) m 

~ar~ -(h)" (5) ~dT (Fe/H) ' (2) 

Here (Fe/H) is the number ratio of Fe to H atoms in the gas, Pp e is the net rate of Fe 
production by SNe II, (H) is the number of H atoms in the gas, and d{B.) in /dt is the rate of 
infall of H atoms into the halo. The time t = corresponds to the big bang. The equation 
governing (H) is: 

rf(H) _ rf(H) as | d(E) iQ 
dt dt dt 
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where d(H)as/rft < is the astration rate. During the infall dominated regime |<i(H) as /<it| <C 
d(H)- m /dt so that (H) ss (H) in . Then for this regime equation (2) can be rewritten as: 



Fe 



^Fe — Ai n Zpe, 



(4) 



where Z Fe = (Fe/H)/(Fe/H) , A Fc = P Fe /[(Fe/H) (H)], and A in = d\n(R)Jdt. Note that 
for Z Fe to increase rapidly, XmZ Fc must be significantly smaller than AF e - 



The general solution to equation (4) for a specific halo is 
Z Fe (t) = exp 









- f \ m (t')dt' 


[ X Fc (t")exp 


f Mt"')dt"' 


- J to 


J to 


J to 



dt", 



(5) 



where to is the onset of astration and Fe production [i.e., Z Fe (t) = for t < to]. The time t 
can be related to the mass of the halo at that time based on the density fluctuation associated 
with the halo (see Fig. 2a and §4). Given the rate of infall A; n and the rate of Fe production 
per H atom in the gas Ap c , Z Fc (t) only depends upon the time t for onset of astration and 
the time t. For the infall phase, we explicitly assume that P Fc /(R) is constant and take 
Afc = (30 Gyr) -1 . This value of Ap e corresponds to the average rate of Fe production by SNe 
II in the Galaxy [~ 1/3 of (Fe/H) over a period of ~ 10 Gyr]. The infall rate A in is not an 
arbitrary function of t but is directly determined from the model of hierarchical structure 
formation discussed in Barkana & Loeb (2001) (see below). The evolution of [Fe/H] in halos 
with infall of big bang debris but without gas expulsion by SNe II can then be calculated in 
a straightforward manner. 

The rate of gas infall governs the growth of the baryonic mass of a halo. We assume 
that a constant fraction Qiy/Q m of the total halo mass is in baryonic matter, where Qb and 
Q m are the fractions of the present critical density contributed by baryonic and all matter, 
respectively. The infall rate is then related to the total halo mass M (mostly in CDM) as 



Ain(t) = 



cHnM dXxvMdz 



dt 



dz dt 



For a flat universe, we have 



dz 
~dl 



-H (i + z )^n m (i + z) 3 + n A , 



(6) 



(7) 



where H = 100/i km s" 1 Mpc^ 1 is the present Hubble parameter and £~2a is the fraction 
of the present critical density contributed by the cosmological constant. We take h = 0.7, 
il b = 0.045, Vl rn = 0.3, and il A = 0.7. Based on the model of hierarchical structure formation 
(see §4), we show M as a function of z for halos associated with la and 2a density fluctuations 
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and 2a halos), respectively, in Figure 2a. Using the function M(z) for a specific halo, 
we calculate A in (t) from equations (6) and (7). Assuming that astration starts (see Fig. 2a 
and §4) when a 1(7 halo reaches M(zq) = 2.8 x 10 8 M Q corresponding to z = 3.8 (to = 1-6 
Gyr) or when a 2a halo reaches M(z ) = 8.4 x 1O 7 M corresponding to z Q = 9.7 (to = 0.49 
Gyr), we numerically integrate equation (4) to calculate the evolution of [Fe/H] = logZ Fc as 
a function of z until z = 1.6 [t = 4.1 Gyr, M(z) = 10 11 M Q ] for the la halo or until z = 4.3 
[t — 1.4 Gyr, M(z) = 10 11 M ] for the 2a halo. The results are shown in Figure 2b. If we 
consider smaller values of M(z ), say M(z ) ~ 1O 5 M , then the onset of astration will be 
moved to higher redshift. The evolution of [Fe/H] for a 2a halo with M(z ) = 7.1 x 10 4 M o 
corresponding to z = 16.8 (t = 0.23 Gyr) is shown as the dot-dashed curve in Figure 2c, 
where the case for M(zq) = 8.4 x 10 7 M Q corresponding to z$ = 9.7 is also shown as the 
solid curve. It can be seen that the dot-dashed curve merges with the solid curve at z ~ 8. 
This is because [Fe/H] quickly approaches a quasi-steady state subsequent to the onset of 
astration (see below). 

Inspection of equation (4) suggests that if A; n (t) ^> \d\ia Z Fe (t)/dt\, then a quasi-steady 
state (QSS) exists such that 

Z Fc (t) » Z^ s (t) = A Fe /A in (t). (8) 
The approach of Z Fe (t) to Z^ ss (t) may be approximated as 



Z Fc (t)^Z^ s (t)-Z^ s (t )exp 



X- m (t')dt' 



t 



(9) 



For t > t + [A in (to)] _1 , the term containing the exponential in equation (9) is insignificant 
and Z Fe (t) approximately equals the QSS value Z^ ss (t) (see Qian & Wasserburg 2003b). 
When we numerically integrate equation (4) to calculate the evolution of [Fe/H] = logZ Fc 
as a function of z, we always use the actual function M(z) for a specific halo to obtain the 
corresponding Aj n (t). For a simple estimate of Aj n (t), we note that the slopes of the M(z) 
curves are \d\n M / dz\ ~ 2.3 within a factor of a few for both la and 2a halos (see Fig. 2a). 
We also note that for z > 1, ft a can be ignored in equation (7) to give 

t w 1 w 17(1 + z)- 3/2 Gyr. (10) 

3 Hoy/fijl + z f 

Using \d\nM/dz\ ~ 2.3 and equations (6) and (10), we obtain 

Ain(f) ~ 10(Gyr/t) 5/3 Gyr" 1 ~ 0.09(1 + z) 5/2 Gyr" 1 (11) 
as a rough estimate for both la and 2a halos. This gives for these halos 

Z£ ss (t) ~ (1/300) (t/Gyr) 5 / 3 ~ 0.4(1 + z)^ 2 , (12) 



-7- 



[Fe/H] QSS = log Z^ s (t) ~ -2.9 - 2.5 log . (13) 

The above estimate for [Fe/H]q SS is shown as the solid curve labeled "quasi-steady state 
estimate" in Figure 2b. 

It can be seen from Figure 2b that for both la and 2a halos, [Fe/H] of the gas initially 
rises very rapidly and approaches the QSS value, which ranges from [Fe/H]Qss ~ —2.8 at 
z = 8 to —1.6 at z — 2. For comparison, the case of evolution from an initial metal-free 
state without infall, Z Fe (t) = \ Fc t (prior to the onset of SNe la), is also shown in Figure 2b. 
This should give an upper bound to Z Fe for all halos and appears to be a reasonable bound 
for the observational data on DLA systems (Wasserburg & Qian 2000b; Qian & Wasserburg 
2003b; see Fig. 5 here). Compared to the case of no infall, the Fe enrichment for la and 2a 
halos is suppressed due to dilution by infall by a factor of ~ A in (i)£ ~ 10(Gyr/t) 2 / 3 , which 
ranges from a factor of ~ 4.5 for t = 3.3 Gyr (z = 2) to a factor of ~ 14 for t = 0.63 Gyr 
(z = 8). 



3. Effects of Infall Cessation 

We now consider that infall of gas ceases at some time for a halo and investigate the 
subsequent chemical evolution of gas in the halo. The metallicity distribution of stars formed 
before and after infall cessation (IC) is studied. The general effects of SNe la on the chemical 
evolution of gas and stars in a halo are also discussed at the end of this section. 



3.1. Chemical Evolution of Gas 

We first consider the chemical evolution of gas in a halo after continuous infall ceases 
sharply at some time tic- The governing equations in this case are still equations (2), (3), 
and (4). But now d(H) [n /dt = 0, Ai n = 0, and the new initial state has ^Fe(^ic) ~ ^R; SS (^ic) 
as determined by the prior evolution with infall. Thus, for t > tic but prior to the onset of 
SNe la (see §3.3), 

Z Fe (t) = Z Fc (t lc ) + \ Fc (t - t lc ). (14) 

Taking tic — 4.1 Gyr (zjc = 1.6) for a la halo and tic — 1-4 Gyr (zic = 4.3) for a 2a halo, 
we show the above evolution of [Fe/H] = log-Zpe for these two halos in Figure 2b. It can 
be seen that the evolution subsequent to infall cessation departs markedly from the QSS 
case. For the 2a halo, the evolution clearly begins to approach the case of no infall. In 
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these calculations, we have assumed that infall ceases sharply at time tic corresponding to 
redshift zic when the mass of a halo reaches M(zic) — 10 11 M Q (see Fig. 2a and §4). If we 
consider larger values of M(z IC ), say M{z\c) = 10 12 M Q , this will only delay the onset of 
uniform growth of Z Fc to lower redshift. The general evolution of [Fe/H] outlined above will 
be unaffected. The evolution of [Fe/H] for a 2a halo with M(zjc) = 10 12 M & corresponding 
to tic = 2.3 Gyr (zic = 2.8) is shown as the dashed curve in Figure 2c. 



3.2. Metallicity Distribution of Stars 

In addition to the effect on the evolution of Z Fe discussed above, the instant cessation of 
infall has another consequence. During the infall phase, the amount of gas in a halo rapidly 
increases as (H) ?a (H)i n for |o?(H) as /<it| <C rf(H); n / 'dt. Subsequent to instant infall cessation, 
(H) steadily decreases due to astration. Thus, (H) is maximal at time tic- As the rate of 
star formation increases with (H), this rate also reaches a maximum at time tic- 

The general history of star formation before and after infall cessation is reflected by the 
metallicity distribution of stars in a halo. This distribution is defined as the number (N) of 
stars formed per unit interval of [Fe/H] = \ogZ Fe , 

dN (lnl0) ^g)^M (15) 



d[Fe/H] v ; (H)* dZ Fe /dt ' 

where (H)* is the average number of H atoms in a star [note that d(R) as /dt < 0]. 

The metallicity distribution of stars formed at t > ti C is simple to obtain. In this case, 
d(R)/dt = ct(H) as /dt and dZ Fc /dt = P Fe /[(Fe/H) (H)] (cf., eqs. [3] and [4]). Assuming that 
the SN II rate is the same as the rate at which massive SN II progenitors are formed (the 
instantaneous recycling approximation), we can relate P Fc to the astration rate d(H) as /<it as 

P Fc = -a(Fe/H) ^, (16) 

where a is a positive dimensionless constant. Thus, for t > tic, 

dZ Fc _ a tt(H) as _ ctln(H) 
dt (H) dt ~ a dt ' 1 ' 

The above equation can be solved to give 

Z Fc — Z Fc (tic) 



(H) = (H) Ic exp 



a 



(18) 
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where (H)ic is the number of H atoms in the gas at t — tic- The metallicity distribution at 
t > tic is then 



dN n ^^Z Fc d(R) (lnlO)(H) IC 



-(lnlO)-^-^ = * - "'"^ exp 



d[Fe/H] v J {E)*dZ Fc a (H) 



Z Fc — Zpcitic) 



a 



(19) 



The above metallicity distribution assumes that only SNe II provide Fe. Note that this 
distribution does not depend on the specific forms of rf(H) as /c?t and P Fe as functions of time 
but only on the ratio of these two quantities. This is true of any closed system as recognized 
by Hartwick (1976) long ago. 

In contrast, the metallicity distribution of stars formed during the infall phase depends 
on the specific forms of d(H) as /<it and P Fc . In general, the evolution of Z Fe and (H) during 
the infall phase must be solved together as explicit functions of time from equations (2) and 
(3) after d(R) as /dt, P Fe , and the infall rate d(H) m /dt are specified for a given halo. We again 
assume P Fe /[(Fe/H) (H)] = A Fe = (30 Gyr) -1 as in §2. Using this assumption and equation 
(16) gives 

fflk = _^( H ). (20) 
dt a y ' y ' 

As a numerical example of estimating a, we assume that SNe II with progenitors of ~ 20- 

25 M produce ~ 0.1 M of Fe as in the case of SN 1987A but no other SNe II produce 

Fe. For a Salpeter initial mass function dN/dM+ oc M~ 2 - 35 over a range of stellar masses 

M* ~ 0.1-100 M Q , a total of ~ 1700 M of stars are formed for each Fe-producing SN II. 

From equation (16) we obtain a ~ 0.1 M /(1O -3 x 1700 M ) ~ 1/17, where we have used an 

Fe mass fraction of ~ 10~ 3 corresponding to the number abundance (Fe/H) . Noting that 

the Fe yields of SNe II, especially the dependence of these yields on the progenitor masses, 

are rather uncertain, we take a = 1/18 as a nominal value and examine the effects of a 

higher value of a — 1/6 later. For a halo with 10 10 M of gas and Af c = (30 Gyr)" 1 , the 

astration rate in equation (20) is 6 and 2M yr _1 for a = 1/18 and 1/6, respectively. These 

astration rates are typically much smaller than the infall rate as assumed in §2. However, to 

be more accurate, we can solve Z Fe and (H) together as functions of time from equations (2), 

(3), and (20). These two functions can then be used to obtain the metallicity distribution 

at t < tic from equations (15) and (20). The metallicity distribution over the whole time 

domain t < tic an d t > tic is shown for a 2a halo with a — 1/18, to — 0.49 Gyr (z = 9.7), 

and instant infall cessation at tic = 1-4 Gyr (zic = 4.3) in Figure 3a (solid curve). It 

can be seen that this distribution has a sharp peak at [Fe/H] = —2.2 and a broad peak 

at [Fe/H] = —1.3. The broad peak occurs in the region t > tic and is easy to understand 

from equation (19). This equation shows that dN/d[Fe/H] has a peak at Z Fe = a, which 

corresponds to [Fe/H] = log a = —1.3 for a = 1/18. 
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The sharp peak in the metallicity distribution of stars occurs at t — tic, the time of 
infall cessation. The origin of this peak can be understood by considering equation (15) in 
the immediate neighborhood of tic (^ic < ^ic < ^ic)- Note that dZ Fc /dt has a discontinuity 
at t = tic but Z Fc and ct(H) as /<it are the same on both sides of this discontinuity. This gives 

(dN/d[Fe/R]) t - c = {dZ Fc /dt) ttc = A Fe 
~ /+ " (dN/d[Fe/R]) ttc (dZ Fc /dt) t - c (dZ Fe /dt) t - c [ } 

As the QSS approximation Z Fc rs Zp^ s holds during the infall phase (t < tic), (dZ-p e / dt) t - ~ 
to zeroth order and R-/+ is singular (see eq. [4]). To first order, 



dZ Fe \ I dZ Fc 



QSS 



dt J .- \ dt 



tic 



d(X Fe /X in ) 
dt 



X Fe ( tic \ 



tic 



2/3 



6 \GyrJ 



(22) 



where we have used the estimate for A; n in equation (11) for la and 2a halos to obtain the 
last result. It follows from equations (21) and (22) that the ratio R-/+ ~ 6(Gyr/ti C ) 2 ^ 3 for 
la and 2a halos is much greater than unity for a wide range of tic- This leads to a sudden 
drop in dN/d[Fe/H] at t = tic an d produces a sharp peak in the metallicity distribution. 
The value of [Fe/H] at the sharp peak depends on t ic . For a 2a halo with t IC = 1.4 Gyr 
(zic = 4.3), this peak is at [Fe/H] ~ log Zpf s (ti C ) 2.2 (see eq. [12] and Fig. 3a). 

To compare the effect of instant infall cessation as assumed above with that of a smooth 
transition for a 2a halo, we have calculated models where d(H) in /dt decays smoothly starting 
at time tic with the functional form exp[— (t — ti C )/r] for t > tic- The results are shown in 
Figure 3a for r = 0.1 and 0.2 Gyr. As can be seen, the main effect of changing r from 
(instant infall cessation) to a few tenths of a Gyr is to increase the fraction of stars formed 
at t > tic- This is because there is more gas available for star formation in the case where 
the infall rate decays smoothly over an extended time. However, the width of the sharp 
peak is increased only slightly from ~ 0.06 dex for r = to ~ 0.14 dex for r = 0.2 Gyr. 
Concerning the total population of stars in the metallicity distribution, we note that only 
10% of all the stars are in the sharp peak for r = and this fraction is even smaller for a 
finite r. For a slowly-decaying infall rate with r ~ 1 Gyr, the sharp peak at t = tic is still 
present but contains only ~ 2% of the stars and becomes a small feature of an otherwise 
smooth metallicity distribution (see Fig. 3b for the case of r = 1 Gyr). 

The bimodal form of the metallicity distribution of stars discussed above applies to all 
halos in general. One critical aspect of this bimodal behavior is the value of [Fe/H] at the 
sharp peak. This value is governed by the time tic for infall cessation, which depends on 
the density fluctuation associated with a halo (see §4). The metallicity distributions for a 
la halo with instant infall cessation (r = 0) and smooth transition to exponential decay of 
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infall (r = 1 Gyr) at a relatively late time of tic = 4.1 Gyr (zic = 1.6) are shown in Figure 
3c [a = 1/18 and t = 1-6 Gyr (z = 3.8) are used]. Compared with Figures 3a and 3b 
for a 2a halo with tic = 1-4 Gyr (zic = 4.3), the sharp peak for the la halo is shifted to a 
higher [Fe/H] of —1.8. To obtain a peak with a lower [Fe/H] of —3 would require tic ~ 0.5 
Gyr corresponding to zic ~ 10 (see eq. [13]). Comparison of Figures 3a and 3b with 3c also 
shows that the sharp peak for the 2a halo contains a much smaller fraction of stars than 
that for the la halo. Consider the case of instant infall cessation as a specific example. For 
the 2a halo, the infall rate at high redshift corresponding to t < t K is large compared to the 
astration rate and only 10% of stars are formed during the infall phase. In contrast, for the 
la halo, the astration rate is no longer small compared to the infall rate at lower redshift 
near the end of the infall phase and 28% of stars are formed during this phase. From the 
model presented here it follows that galaxies in general should have a bimodal distribution 
of stars with a relatively small number of stars defining a sharp peak at low [Fe/H]. This 
peak reflects the state when gas infall ceased or greatly decreased. Note that we have not 
addressed the location of stars in a galaxy when discussing their metallicity distribution. 
The bimodal metallicity distribution is more likely to describe stars in the halo than those 
in the disk of a galaxy. 

As discussed above, the metallicity distribution of stars in a halo is calculated by speci- 
fying AFe, a, to, he, and the form of infall cessation. The infall rate Ai n used in the calculation 
for t < tic is determined from the model of hierarchical struture formation as in §2. The 
effects of the form of infall cessation (instant or exponential decay) on the metallicity distri- 
bution have been discussed above. To examine the effects of a and tic, we have calculated 
the metallicity distribution for a 2a halo with a = 1/6, to = 0.49 Gyr (zq = 9.7), and instant 
infall cessation at tic = 2.3 Gyr (zic = 2.8). This is shown as the solid curve in Figure 
3d. Compared to the solid curve in Figure 3a [a = 1/18 and tic = 1-4 Gyr (zic = 4.3)], 
the sharp peak in Figure 3d is shifted to [Fe/H] = —1.9 ~ log Z^, ss (tic) due to a later tic 
and the broad peak (solid curve) in Figure 3d is shifted to [Fe/H] = logo; = —0.8 due to 
a larger a. In addition, 15% of all the stars represented by the solid curve in Figure 3d 
have —0.5 < [Fe/H] < while the corresponding fraction is only 0.3% for the solid curve in 
Figure 3a. For completeness, we note that t for the onset of astration has little effect on 
the metallicity distribution as most of the stars are formed close to t = t K for t < tic- 

It is evident from inspection of Figure 3 that if a protogalactic halo were to expel 
its gas immediately following the instant cessation of infall, an assemblage of stars with 
nearly identical [Fe/H] would be left behind. This assemblage of stars may be accreted 
by a larger system and become a globular cluster of the larger system. In this model, the 
sharp distribution of [Fe/H] for the stars in a globular cluster would not be the result of a 
rapid burst of star formation but rather the result of the following two factors. First, for a 
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protogalactic halo with ongoing infall and astration over an extended period, the rate of star 
formation increases rapidly with time and reaches a maximum at the time tic f° r instant 
infall cessation when there is the largest amount of gas in the halo. This causes nearly all 
the stars formed during the infall phase to have the metallicity of the gas at t — tic- Second, 
the rapid loss of gas subsequent to infall cessation terminates further astration. Our model 
relates the value of [Fe/H] for the globular cluster to t — tic, and therefore, can be used to 
estimate the formation time (relative to the big bang) of the cluster based on its [Fe/H] if 
the initial metallicity of the infalling gas is negligible. 



3.3. General Effects of SNe la 

To illustrate the general effects of SNe la, we assume that these supernovae start to 
contribute Fe when [Fe/H] = — 1 is reached at time t'. The rate of Fe production by SNe 
la is taken to be twice that for SNe II so that a solar Fe abundance can be achieved over 
a period of ~ 10 Gyr as in the Galaxy. The total Fe production rate per H atom in the 
gas then increases from Af c at t < t" 1 " to 3A Fc at t > tf. Consider halos with instant infall 
cessation occurring at t — tic prior to the onset of SNe la (i.e., > tic). The evolution of 
Z-p c at t > for these halos is given by 

Z Fc (t) = Z Fe (tt) + 3A Fc (t - t f ), (23) 

where Z-p e {p) = 0.1 ([Fe/H] = —1) is taken to correspond to the onset of SNe la. The 
above evolution of [Fe/H] = \ogZ Fc is shown for a la halo with t = 1.6 Gyr (z = 3.8) and 
tic = 4.1 Gyr (z ic = 1.6) and for a 2a halo with t = 0.49 Gyr (z = 9.7) and t IC = 1.4 Gyr 
(zic = 4.3) in Figure 2b. The evolution of [Fe/H] at t > v for the case of no infall is also 
given by equation (23) and shown in Figure 2b. The only effect of SNe la on the evolution 
of [Fe/H] for all halos is to increase the rate of Fe enrichment at t > tf. 

As in the case of infall cessation, the sudden increase of the Fe production rate per H 
atom in the gas at the onset of SNe la leads to a sharp drop of dN/d[Fe/H] at [Fe/H] = — 1. 
This is shown for a 2a halo with a = 1/18 and smooth transition to exponential decay of 
infall (r = 1 Gyr) at tic = 1.4 Gyr (z IC = 4.3) in Figure 3b and for another 2a halo with 
a — 1/6 and instant infall cessation at tic = 2.3 Gyr (z IC = 2.8) in Figure 3d. For the 
cases where only SNe II provide Fe, the metallicity distribution of stars has a broad peak 
at [Fe/H] = logo; = —1.3 for Figure 3b and —0.8 for Figure 3d (solid curves). If we include 
SN la contributions for [Fe/H] > — 1, this broad peak is truncated and a new broad peak is 
formed at [Fe/H] = log(3ai) = —0.8 for Figure 3b and —0.3 for Figure 3d (dashed curves). 
This change to the metallicity distribution is caused by the more rapid increase of [Fe/H] 
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with additional contributions from SNe la. The above effect of SNe la on the metallicity 
distribution of stars should apply to all halos. For completeness, we note that the final halo 
mass for the 2a halo shown in Figure 3b is 5.5 x 10 11 M corresponding to a baryonic mass of 
8.3 x 10 10 Mq. In this case 89% of the baryonic matter has been consumed in star formation 
prior to the onset of SNe la. In comparison, the final halo mass for the 2a halo shown in 
Figure 3d is 10 12 M corresponding to a baryonic mass of 1.5 x 10 11 M . In this case 45% 
of the baryonic matter has been consumed in star formation prior to the onset of SNe la. 



In the formal discussion of §2 and §3 on the chemical evolution of a halo, we used the 
times for the onset of astration (to) and infall cessation (tic) and the infall rate (Am), all of 
which are related to the growth of the halo as determined by the CDM model of hierarchical 
structure formation. We now explicitly explore how this cosmological input is determined. 
The physical conditions governing the general evolution of a halo are closely related to the 
mass (mostly in CDM) of the halo. For a halo associated with an na density fluctuation (an 
na halo), its mass M at redshift z > 1 can be estimated from 



where a(M) is the standard deviation of the present power spectrum of overdensity on mass 
scale M. Using the version of equation (24) presented in the Appendix and taking a(M) 
from Figure 5 of Barkana & Loeb (2001), we show M as a function of z for la and 2a halos, 
respectively, in Figure 2a. The function M(z) for a specific halo was used to calculate the 
infall rate A m from equations (6) and (7) in §2. In addition, when discussing the effects of 
infall cessation in §3, we assumed that infall ceases at time tic corresponding to redshift zic 
at which the mass of a halo reaches M(zi C ) = M IC . For Mi C = 1O 11 M , z IC = 1.6 or 4.3 
for a la or 2a halo, respectively (see Figure 2a). If Mj C = 10 12 M is assumed, then infall 
cessation is delayed to zjc = 0.8 or 2.8 for a la or 2a halo, respectively. The redshifts at 
which infall cessation occurs for other na halos are given for Mjc = 10 11 M and 1O 12 M , 
respectively, by the corresponding curves in Figure 4. The range of M\q = 10 11 to 10 12 M 
is in accord with the halo masses of typical galaxies. 

A critical physical parameter for the onset of astration in a halo is the virial temperature 
T vir of its gas. For z > 1 this temperature is related to the halo mass M as 



4. Cosmological Input for Halo Evolution 



na(M) w 1.33(1 + z), 



(24) 




(25) 
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where \x is the mean molecular weight and fi = 1.22 or 0.6 for a neutral or ionized gas of 
primordial composition, respectively. Following Barkana & Loeb (2001), we assume that 
astration starts in a halo when the gas reaches a threshold virial temperature T vil . . At a 
given redshift z, the mass M = M that is required for a halo to reach T vir = T vil . can be 
obtained from equation (25). Using this equation (see Appendix), we show Mq as a function 
of z for T viri0 = 300 K (// = 1.22) and T viri0 = 10 4 K (// = 0.6), respectively, in Figure 2a. 
The redshift z = z at which the mass M(z) of an na halo matches M (z) then corresponds 
to the onset of astration (t ) in this halo. For T virj0 = 10 4 K, z = 3.8 or 9.7 for a la or 2<r 
halo, respectively (see Figure 2a). If T vil:0 = 300 K is assumed, then the onset of astration 
occurs earlier at z = 7.3 or 16.8 for a la or 2a halo, respectively. The redshifts at which 
astration starts for other na halos are given for T vir)0 = 300 or 10 4 K, respectively, by the 
corresponding curves in Figure 4. 

The value of T vir)0 for the onset of astration depends crucially on the available coolant. 
In the presence of cooling by H2 molecules, T v ; r o = 300 K is appropriate. If H2 molecules 
are dissociated (e.g., Couchman & Rees 1986; Haiman, Rees, & Loeb 1997; Ciardi, Ferrara, 
& Abel 2000), then T v ; r o = 10 4 K is required for efficient cooling by atomic species. While 
the choice of T vil . = 300 or 10 4 K does not affect the long-term chemical evolution of a halo 
(see Fig. 2c), it can make an important difference in the early evolution of the halo. This 
difference concerns the possibility for a single SN II to expel all the gas, thereby disrupting 
the chemical evolution of a halo. This possibility can be assessed by comparing the explosion 
energy of ~ 10 51 erg for an SN II with the binding energy -Ebi, g as of the gas in the halo. For 
z > 1 this binding energy is related to the halo mass M as 

£bi Kas « 4.31 x 10 52 ( f ) ' ( erg. (26) 

We consider a simple scenario where all the gas is expelled when -Ebi, ga s f° r a halo is less 
than 10 51 erg. Using equation (26) (see Appendix), we show the halo mass M = for 
Ebi,gas — 10 51 erg as a function of z in Figure 2a. The redshifts at which different na halos 
reach -Ebi.gas = 10 51 erg are given by the corresponding curve in Figure 4. As can be seen 
from this figure, if T vir = 300 K is assumed for the onset of astration, then the chemical 
evolution of all halos will be disrupted by SNe II until -Ebi.gas = 10 51 erg is reached. In 
contrast, for T vilj0 = 10 4 K, such disruption cannot occur as -Ebi.gas = 10 51 erg is always 
reached prior to the onset of astration. 
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5. Comparison with Data on DLA Systems 

We now compare our chemical evolution model with the data on [Fe/H] for DLA systems. 
Using T v i rj o = 10 4 K for the onset of astration, M\c = 10 11 M for instant infall cessation, 
and [Fe/H] = —1 for the onset of SNe la, we show the evolution of [Fe/H] as a function 
of z for 0.67(7, la, and 2a halos along with the data on DLA systems (Prochaska et al. 
2003) in Figure 5 (see also Qian & Wasserburg 2003b). The case of no infall is also shown 
in this figure. It can be seen that the chemical evolution model presented here provides 
a good description of the data. For any given z the model gives a wide range of [Fe/H] 
corresponding to three different evolutionary stages of DLA systems: sharp initial rise after 
the onset of astration, quasi-steady state due to competition between Fe production and 
dilution by infall, and uniform growth subsequent to infall cessation. The evolutionary 
stage of a DLA system at a given z depends on the density fluctuation associated with the 
system. The large dispersion in [Fe/H] for DLA systems at any given z is then accounted for 
by the different density fluctuations associated with these systems. A detailed quantitative 
assessment of the model would require many more DLA systems with accurate measurements 
of [Fe/H] at a given z (see Qian & Wasserburg 2003b for an illustrative example). The model 
also provides an explanation for the baseline enrichment of [Fe/H] ~ —3 observed for DLA 
systems. Subsequent to the onset of astration, it takes only ~ 30 Myr to reach [Fe/H] ~ —3 
for Ap c = (30 Gyr)" 1 . Due to this rapid initial rise of [Fe/H], the probability to find DLA 
systems with [Fe/H] below ~ —3 is very low. The baseline enrichment of [Fe/H] ~ —3 for 
DLA systems is similar to the abundance level in the IGM that was attributed to a prompt 
inventory of metals suggested by other arguments (e.g., Wasserburg & Qian 2000a; Qian 
& Wasserburg 2002; Qian, Sargent, & Wasserburg 2002). Based on the model presented 
here, this postulated prompt inventory is not required to explain the data on DLA systems 
although it is also not in conflict with these data. 



6. Discussion and Conclusions 

We presented a model for the chemical evolution of the gas and stars in halos as a 
function of redshift. This model assumes that a normal population of stars including SN II 
progenitors can form from big bang debris and that SNe II will provide metals such as the 
elements from the Fe group down to O. If there is a prompt inventory of metals corresponding 
to [Fe/H] < —3 in the IGM, this will simply provide a baseline enrichment for the initial state 
and will not alter the results. The infall of gas into a halo is assumed to closely follow that 
of CDM. The infall rate and the general evolution of the halo are prescribed by the model of 
hierarchical structure formation based on the density fluctuation associated with the halo. 
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For systems that are not disrupted by SN II explosions, the chemical evolution of the gas and 
stars is explicitly determined by the infall rate as compared with the astration rate and the 
corresponding rate of metal production by SNe II per H atom in the gas. Once the astration 
rate and the metal production rate are specified, the only "free" parameters in this model 
are the times (or states) at which astration starts and infall ceases for a halo. For all systems 
not disrupted by SN II explosions, if there is a cessation of gas infall, the metallicities of 
stars follow a bimodal distribution. This distribution is characterized by a sharp peak at 
the value of [Fe/H] corresponding to the time of infall cessation and by a broad peak at a 
higher value of [Fe/H] corresponding to the subsequent period of astration during which the 
bulk of the remaining gas forms stars. The fraction of stars in the sharp peak depends on 
how rapidly the infall cessation takes place. While it is not clear how observational bias may 
affect the proper metallicity distribution of stars, we note that the distribution observed 
for the Galactic halo stars (e.g., Fig. 5 in Christlieb 2003) and the bimodal distribution 
presented here show a strong resemblance. Based on our model, a distinct narrow peak at 
low [Fe/H] for the Galactic halo stars would signify the time when infall started to decrease 
for the Galaxy. We also discussed the general effects of SNe la on the chemical evolution of 
the gas and stars in halos. For [Fe/H] > — 1 the production rate for the Fe group elements 
is increased due to additional contributions from SNe la. Consequently, the broad peak in 
the metallicity distribution of stars corresponding to the case where only SNe II provide Fe 
is truncated at [Fe/H] ~ —1 and a new broad peak is formed at a higher [Fe/H] value. 

We also considered the physical conditions governing the evolution of a halo. The mass 
M of a halo is an important parameter and its growth as a function of redshift is determined 
by the model of hierarchical structure formation based on the na density fluctuation asso- 
ciated with the halo. We assumed that the cessation of infall for both baryonic and dark 
matter corresponds to the redshift at which the halo mass reaches M\q = 10 11 to 1O 12 M . 
This range of M\q is in accord with halo masses of typical galaxies. The other two important 
parameters of a halo are the virial temperature T vir and the binding energy E^ gas of the 
gas. We assumed that astration starts in a halo when a critical virial temperature T viri0 
corresponding to a halo mass M is reached. In the presence of cooling by H 2 molecules, 
T v irfl = 300 K is appropriate and M ranges from ~ 10 5 to 10 6 M for different na halos. If 
H 2 molecules are dissociated, then T vir = 10 4 K is required and M ranges from ~ 10 8 to 
10 9 Mq. While the assumed value of T v i r o does not affect the long-term chemical evolution of 
a halo, it makes an important difference concerning the effect of SN II explosions on the early 
evolution of the halo. For T v ; ri o = 10 4 K, astration starts in a halo after the halo mass reaches 
M bi corresponding to a gas binding energy of -Ebi,gas = 10 51 erg. In this case, the chemical 
evolution of the gas and stars in the halo will not be disrupted by SN II explosions and is 
summarized in the preceding paragraph. In contrast, for T viri0 = 300 K, astration starts in a 
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halo before the halo mass reaches M hi and all low-mass halos with M < M < M hi will have 
their gas dispersed into the IGM by a single SN II. If such low-mass halos indeed form, then 
their gas expulsion by SNe II may significantly enrich the IGM at very high redshift. For 
example, at z — 10, halos with M < M < My are associated with 1.3-1. 7a density fluctua- 
tions (see Fig. 4) and a fraction F{M < M < M bi ) = y/{2/ir) exp(-x 2 /2)dx « 10% of 
all baryonic matter is being processed through such halos. Taking Z Fc ~ Z^ s ~ 10~ 3 (see 
eq. [12]) for the gas expelled from such halos, we estimate that their gas expulsion would 
result in an average metallicity of Z^ M ~ F(M < M < M bi )zgf s ~ 10~ 4 corresponding to 
[Fe/H]i G M ~ —4 in the IGM at z — 10. We expect that T virj0 for the onset of astration would 
increase from 300 to 10 4 K as various radiation sources turned on during the evolution of 
the universe. So the range of low-mass halos that would be disrupted by SN II explosions 
would decrease with time and eventually vanish. The question of whether such low-mass 
halos may actually provide significant metal enrichment of the IGM at very high redshift 
will be discussed in a subsequent paper that considers aspects regarding dissociation of H 2 
molecules. 

We proposed a possible mechanism of globular cluster formation based on the sharp 
peak in the metallicity distribution of stars. This peak corresponds to the time of infall 
cessation. If infall ceases rather suddenly for a protogalactic halo and its gas is rapidly 
expelled afterwards, an assemblage of stars with nearly identical [Fe/H] will be left behind. 
This assemblage of stars may be accreted by a larger system and become a globular cluster 
of the larger system. To estimate the masses and metallicities of such globular clusters, we 
assume T vir = 10 4 K for the onset of astration. In this case the gas in a halo will be enriched 
by SNe II without disruption subsequent to the onset of astration. During the infall phase 
the amount of infalling gas increases rapidly on a timescale of Ar 1 ~ 0.04[10/(1 + z)f/ 2 Gyr 
(see eq. [11]). For an astration rate of <i(H) as /<it = — (AF e /a;)(H) with Af c = (30 Gyr) -1 and 
a = 1/18, the fraction of the halo mass M in stars at redshift z is /* and can be estimated as 
/* ~ (fi 6 /fi m )(A Fc /a)A in 1 ~ 4x 10~ 3 [10/(l+2:)] 5/2 . If infall ceases at redshift z, then the sharp 
peak in the metallicity distribution contains a total mass ~ f+M of stars with essentially 
the same [Fe/H] as given approximately by [Fe/H]q SS for this z (see eq. [13]). The assumed 
value of T virj0 = 10 4 K corresponds to M = 8.1 x 10 7 and 2.6 x 10 8 M Q for z = 10 and 4, 
respectively (see Fig. 2a). So in the above scenario, the globular clusters would have masses 
exceeding /*Mo ~ 3 x 1O 5 M and [Fe/H] ~ —3 if they were formed at z — 10 while they 
would have masses exceeding /*M ~6x 10 6 M and [Fe/H] ~ —2.1 if they were formed at 
z — 4. No small globular clusters with masses of ~ 10 4 M Q are formed in this scenario unless 
the larger aggregates are fragmented. We note that the masses derived for Galactic globular 
clusters are rather uncertain and range from ~ 300 to ~ 10 7 M Q with a mean of ~ 2 x 10 5 M Q 
(Mandushev, Spassova, & Staneva 1991; Richer et al. 1991). The metallicity distribution 
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of Galactic globular clusters is bimodal with peaks at [Fe/H] ~ —1.6 and —0.5, respectively 
(Armandroff & Zinn 1988). Approximately 10% of the clusters have [Fe/H] < —2. The 
masses and metallicities of globular clusters in our model are broadly consistent with those 
of the Galactic population. A detailed comparison requires considerations of possible ranges 
for Af C! Ai n , and a. We also note that stars in an individual Galactic globular cluster typically 
have a very narrow range of [Fe/H], which can be explained by our model. However, there 
is a considerable range in [Fe/H] for the most massive Galactic globular cluster oj Centauri 
(Suntzeff & Kraft 1996; Norris, Freeman, & Mighell 1996) with a mass of 7 x 10 6 M (Richer 
et al. 1991). The range of [Fe/H] for uo Centauri may have resulted from further astration 
subsequent to infall cessation based on our model. As an SN II cannot unbind the gas in 
the protogalactic halo associated with such a large cluster and some additional mechanism 
for gas loss must be involved, a more extensive evolution of [Fe/H] would be expected. The 
gas loss may have to occur by e.g., passing through the Galactic disk. 

If we assume T v i r o = 300 K for the onset of astration, then a different scenario of globular 
cluster formation may occur. This scenario is possible if H 2 molecules are not dissociated and 
the IGM is not reionized so we here only consider it for very high redshift, say z ~ 10. For 
z = 10, M = 1.5 x 10 5 M for the assumed value of T virfi = 300 K and M bi = 9.9 x 10 6 M 
(see Fig. 2a). A low-mass halo with M < M < Mu has started astration but its gas will be 
expelled by a single SN II. As the metals produced by the SN II are expelled along with the 
gas, any metals in the stars already formed cannot be due to this SN II but must have come 
from the infalling gas. In other words, the metallicities of these stars represent the metal 
inventory of [Fe/H] < —3 provided to the IGM at z ~ 10 by sources at even higher redshift. 
The assemblage of stars after the SN II explosion resembles a globular cluster with a mass 
of ~ f*M ~ 500 M to ~ /*Mbi ~ 3 x 1O 4 M . We note that Galactic globular clusters in 
this mass range typically have [Fe/H] much higher than —3. Whether globular clusters in 
the same mass range but with [Fe/H] < —3 exist remains to be seen. The possibility that 
low-mass halos might form at z < 10 from an enriched IGM also remains to be explored. 

We would like to acknowledge the anonymous referee for one of our earlier studies on 
cosmochemical evolution (Qian & Wasserburg 2003a) for being supportive and for suggesting 
that we might pay attention to infall. Good advice sometimes does fall in. We also want to 
thank the anonymous referee for the present work and the Scientific Editor, Brad Gibson, for 
helpful suggestions. This work was supported in part by DOE grants DE-FG02-87ER40328, 
DE-FG02-00ER41149 (Y. Z. Q.) and DE-FG03-88ER13851 (G. J. W.), Caltech Division 
Contribution 8909(1111). 
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A. Equations for Halo Evolution 

In this appendix, we summarize the equations that are used to calculate the mass M, 
the virial temperature T v ; r , and the gas binding energy i?bi,gas of a halo for any redshift z. 
While these equations were discussed in many works, our discussion here closely follows that 
of Barkana & Loeb (2001). We consider a flat universe with Q m = 0.3, Vl b = 0.045, f2 A = 0.7, 
and h = 0.7. 

The evolution of the mass M for an no halo is prescribed by 

1.686^| = na(M), (Al) 

where D(z) is the growth factor (Peebles 1980) and a(M) is the standard deviation of 
the present power spectrum of overdensity on mass scale M. The growth factor can be 
approximated as (e.g., Eisenstein & Hu 1998) 

U ra 4/7 -^+[l + (^/2)][l + (fii/70)]' 1 } 

where 

n* = <Ui + *) 3 (A3 ) 

fi m (i + ^) 3 + fi A ' 1 ; 

Q1 = n m (i + z)* + n A (A4) 

Note that the proportionality relation in equation (A2) is sufficient to determine the ratio 
D(0)/D(z) in equation (Al). For z > 1, equation (Al) can be simplified to give equation 
(24) for the adopted cosmological parameters. 

The virial temperature T vir of the gas in a halo of mass M at redshift z is 

where A c = 18tt 2 - 82fi A - 39(fi A ) 2 ( cf - Br yan & Norman 1998). The binding energy £ bijgas 
of the gas is 

J? KAK m53^ ( M \ 5/3 fQ m A c V /3 fl + z\ f . R . 

Ebi sa s = 5.45 x 10 —— — — — — h erg. (A6) 

,g Q m \10 8 h- l M & ) \n z m 187T 2 J \ 10 J B K J 

For z > 1, equations (A5) and (A6) can be simplified to give equations (25) and (26), 
respectively, for the adopted cosmological parameters. 
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Fig. I. — Schematic diagram of the closed-system model for (Fe/H) in DLA systems. The 
baseline (Fe/H)p is the value corresponding to a postulated prompt inventory in the IGM. 
Normal astration and Fe production by SNe II start in a DLA system at time t* after the big 
bang (t — 0). Different t* values for different DLA systems result in a wide range of (Fe/H) 
at a fixed z. For t* = 0, there is a maximum increase in (Fe/H) (trajectory A). 

Fig. 2. — (a) Evolution of the halo mass M as a function of z calculated for la (dot-dashed 
curve) and 2a (long-dashed curve) halos using the model of hierarchical structure formation 
discussed in Barkana & Loeb (2001). For a specific halo, this evolution determines the infall 
rate. A simple estimate of the infall rate for la and 2a halos can be obtained by taking 
a typical slope of \d\a.M/dz\ ~ 2.3. It is assumed that infall ceases when M reaches Mi©. 
For Mic = 10 11 M Q (horizontal solid line labeled "Mi© = 10 11 M©"), this occurs at redshifts 
zjc = 1.6 and 4.3 for la and 2a halos, respectively, while for M IC = 10 12 M Q (horizontal 
solid line labeled "M IC = 10 12 M Q "), this occurs at redshifts zic = 0.8 and 2.8 for la and 
2a halos, respectively. The time at which the virial temperature of the gas reaches a critical 
value T vir)0 is taken as the onset of astration in a halo. The halo mass M required to reach 
T vir is shown as a function of z for T vir = 300 K [solid curve labeled "M (T vir o = 300 K)"] 
and T vir! o = 10 4 K [solid curve labeled "M (T vil . j0 = 10 4 K)"], respectively. The redshift Zq 
at which the curve M (z) intersects the curve M(z) for a halo corresponds to the onset of 
astration. For T vir)0 = 300 K, z = 7.3 and 16.8 for la and 2a halos, respectively, while for 
T V irfl = 10 4 K, z = 3.8 and 9.7 for la and 2a halos, respectively. It is assumed that a single 
SN II can expel all the gas from a halo before the binding energy of the gas -Ebi.gas reaches the 
explosion energy of 10 51 erg for an SN II. The halo mass M bi required to reach E h ^ gas = 10 51 
erg is shown as a function of z (short-dashed curve labeled "M bi "). (b) Evolution of [Fe/H] 
for la (dot-dashed curve) and 2a (long-dashed curve) halos as functions of z. For both 
halos, it is assumed that the onset of astration corresponds to T v ; ri o = 10 4 K, instant infall 
cessation occurs for Mi© — 10 11 M©, and SNe la start to contribute Fe at [Fe/H] = — 1. The 
solid curve labeled "no infall" is for a system with no infall. For all cases, the Fe production 
rate per H atom of the gas is taken to be Ap e = (30 Gyr) -1 when only SNe II provide Fe 
([Fe/H] < -1) and 3A Fe when both SNe II and la provide Fe ([Fe/H] > -1). Using the simple 
estimate of the infall rate, the case where the Fe production by SNe II only is balanced by the 
dilution due to infall is shown as the solid curve labeled "quasi-steady state estimate." The 
overall evolution of [Fe/H] for a 2a halo is indicated by arrows. Note the sharp initial rise of 
the evolutionary trajectories for both la and 2a halos subsequent to the onset of astration. 
These trajectories then lie close to the quasi-steady state estimate until infall cessation (z\c)- 
They rapidly approach the case of no infall after infall cessation (z < z\c). (c) Effects of 
T vir o and M IC on the evolution of [Fe/H] for a 2a halo. The solid curve is the case shown 
in (b). The dot-dashed curve assumes a different T vil . )0 of 300 K and the short-dashed curve 



-23- 



assumes a different M\ C of 10 12 M Q . Note that these three cases share a common period of 
quasi-steady state evolution. 

Fig. 3. — (a) Metallicity distribution (in arbitrary units) of stars in a 2a halo. It is assumed 
that astration starts at z = 9.7 corresponding to T virj0 = 10 4 K and infall starts to decrease 
exponentially on a timescale r at zjq = 4.3 corresponding to M IC = 10 11 M Q . Only SNe 
II provide Fe at a rate of Af c = (30 Gyr) -1 per H atom in the gas. The parameter a for 
the astration rate is taken to be 1/18. The sharp peak in the metallicity distribution at 
[Fe/H] = —2.2 represents the stars formed during the infall phase (z > zic) while the broad 
peak at [Fe/H] = logo; = —1.3 represents those formed subsequently (z < zic) until the gas 
is exhausted by astration. The rapid drop of the metallicity distribution at z — zic is due 
to the sharp change in dZp e /dt between the regimes with and without infall. The fraction 
of stars in the sharp peak is 10% for the case of instant infall cessation (r = 0, solid curve). 
This fraction is even smaller for finite r [dashed (r = 0.1 Gyr) and dot-dashed (r = 0.2 
Gyr) curves] as more stars are formed after infall starts to decrease, (b) Solid curve same as 
(a) but for r = 1 Gyr. The sharp peak at [Fe/H] = —2.2 is still present but contains only 
2% of the stars and is subdued as compared with the cases of rapid infall cessation (r < 1 
Gyr). The dashed curve assumes a total Fe production rate of 3Ap e per H atom in the gas 
for [Fe/H] > — 1 due to additional Fe contributions from SNe la. The sharp drop of the 
metallicity distribution at [Fe/H] = — 1 corresponds to the onset of SNe la. Note that a new 
broad peak is formed at [Fe/H] = log(3a) = —0.8. Prior to the onset of SNe la 89% of the 
baryonic matter has been consumed in star formation, (c) Same as (a) but for a la halo. 
Astration starts at z = 3.8 corresponding to T vir;0 = 10 4 K and infall starts to decrease 
exponentially on a timescale r (solid curve: r = 0, dashed curve: r = 1 Gyr) at z\ G = 1.6 
corresponding to M\q = 1O 11 M . Note that the sharp peak corresponding to the time of 
infall cessation is shifted to [Fe/H] = —1.8 as compared to [Fe/H] = —2.2 for the 2a halo, 
(d) Same as (a) and (b) but for zic = 2.8 corresponding to M\c = 1O 12 M , r = 0, and 
a = 1/6. Note that the sharp peak corresponding to the time of infall cessation is shifted 
to [Fe/H] = —1.9 as compared to [Fe/H] = —2.2 in (a) and (b) while the broad peak of the 
solid curve is shifted to [Fe/H] = —0.8 as compared to [Fe/H] = —1.3 in (a) and (b). Also 
note that the new broad peak (dashed curve) due to additional Fe contributions from SNe 
la is shifted to [Fe/H] = —0.3 as compared to [Fe/H] = —0.8 in (b). Prior to the onset of 
SNe la 45% of the baryonic matter has been consumed in star formation. 

Fig. 4. — Stages of halo evolution. The redshift at which astration starts in an na halo is 
given for T vir = 300 and 10 4 K, respectively, by the correspondingly labeled solid curves. 
The redshift at which an na halo reaches £t>i,gas = 10 51 erg is given by the correspondingly 
labeled dot-dashed curve. The redshift at which infall ceases for an na halo is given for 
Mic = 10 11 and 10 12 M , respectively, by the correspondingly labeled dashed curves. Note 
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that if the onset of astration corresponds to T viri0 = 300 K, then a single SN II can expel all 
the gas from a halo, thereby disrupting its evolution, until -Ebi.gas = 10 51 erg is reached. In 
contrast, if the onset of astration corresponds to T vil . = 10 4 K, then such disruption cannot 
occur as -Ebi.gas — 10 51 erg is always reached prior to the onset of astration. It is expected 
that T v i r o would increase from 300 to 10 4 K as H2 molecules were dissociated during the 
evolution of the universe. Thus, the range of halos that are disrupted by SN II explosions 
would decrease with time and eventually vanish as indicated by the curvy arrow. 

Fig. 5. — Evolution of [Fe/H] as functions of z for 0.67a (short-dashed curve), la (dot-dashed 
curve), and 2a (long-dashed curve) halos compared with the data (squares) for 96 DLA 
systems (Prochaska et al. 2003). The evolution of [Fe/H] for all three halos is calculated 
using the same assumptions as in Figure 2b and the dot-dashed, long-dashed, and solid 
curves are the same as in that figure. The prompt inventory proposed earlier (horizontal 
long-dashed line) is shown for reference. 
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